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ABSTRACT 

We consider the minimum mass planet, as a function of radius, that is capable of opening a 
gap in an a-accretion disc. We estimate that a half Jupiter mass planet can open a gap in a 
disc with accretion rate M < lO _8 M /yr for viscosity parameter a — 0.01, and Solar mass 
and luminosity. The minimum mass is approximately proportional to M O A8 a o s M^ A2 L~ om . 
This estimate can be used to rule out the presence of massive planets in gapless accretion discs. 
We identify two radii at which an inwardly migrating planet may become able to open a gap 
and so slow its migration; the radius at which the heating from viscous dissipation is similar to 
that from stellar radiation in a flared disk, and the radius at which the disc becomes optically 
thin in a self-shadowed disc. In the inner portions of the disc, we find that the minimum 
planet mass required to open a gap is only weakly dependent on radius. If a migrating planet 
is unable to open a gap by the time it reaches either of the transition radii, then it is likely to 
be lost onto the star. If a gap opening planet cuts off disc accretion allowing the formation of 
a central hole or clearing in the disc then we would estimate that the clearing radius would 
approximately be proportional to the stellar mass. 
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1 INTRODUCTION 

Recent observations have identified young (1-3 Myr old) stars that 
have inner clearings in t he dust distribution as inferred from their 
IRS spectra (C oKuTau/4.lD'Alessio et al]|2005l ; GM Aur, TW Hya 
and DM Tau , ICalvet et al.ll200l l2005 | rand brown dwarf candi- 
dates L316 and L30003 in IC 348, IMuzerolle et all 120061) . These 
disks have been dubbed "transitional disks" as they represent the 
stage in which circumstellar disks are probably disappearing and 
in which massive planets are thought to be forming. In some cases 
there is still accretion on to the central star; DM Tau and GM Aur 
at rates of 2 x 10~ 9 and lO _8 M0/yr, respectively. In other cases 
there is no evidence for accretion, i.e., CoKuTau/4, even though 
its outer disk resembles that of other T-Tauri stars. Multi-band 
IRAC/MIPS photometric measurement from imaging surveys have 
been used to identify dozen s of other transition disk candidates (e.g. 
ISicilia-Aguilar et all2006h . Candidates are then observed with IRS 
allowing accurate meas urements o f edge wall heights, and temper- 
atures and so radii (e.g. lKimll2007h . 

Two dominant approaches exist toward predicting and ac- 
counting for holes or clearings in young circumstellar disks. The 
first is planet formation followed by an op ening of a ga p and 
subsequent clearing of an inner disk (e.g. lOuillen et al.l 



ip and 
120041 : 



* rge21@pas.rochester.edu 
f aquillen@pas.rochester.edu 
| jaehong@pas.rochester.edu 



IVarniereetaLI 120061 : ICrida and Morbidellil |2007|). Alternatively , 
photo-evaporation (e.g. Clarke et alj|200lt I Alexander et alj|2006l) 
may clear the disc from the inside out. A third possibility is that 
smaller dust particles ar e preferentia lly destroyed at a particular ra- 
dius due to turbulence (Ciesla 2007), within the context of a dead 
zone model. Photo-evaporatio n models ca n account for large disk 
clearings near luminous stars dK~imll2007l) and non-accreting sys- 
tems such as CoKuTau/4. However they are not sufficiently so- 
phisticated that they can account for clearings in disks hosted by 
low lu minosity brown dwarfs (UV flux is too low. IMuzerolle et al.l 
2006)), disks with wide radial gaps and inner disks (such as GM 
Aur) or disks that are truncated in the dus t distribution but continue 
to acc rete (such as DM Tau). Recently, IChiang and Murray- Ctevl 
d2007l) have suggested that stellar X-rays would be able to ionise 
the inner wall of a gap enough to trigger accretion via the magneto- 
rotational instability. The disc would then clear from the inside out, 
at an accelerating rate. Models involving a planet are likely to be 
more versatile but because of the complexity of planet/disk interac- 
tions this scenario has been little explored. In this paper we focus 
on the possibility that massive planets open gaps in disks and ex- 
plore what mass planets are capable of opening gaps in different 
accretion disks. 

Because a planet drives density waves into a disk 
dLin an d Papaloizou 1 19791 ; [Goldreich and Tremainel Il980l) if the 
planet is sufficiently massive it can o vercome the effect of vis- 
cosity and open a gap in the disk |L in and Papal oizou! 1 19931 : 
IWard and Hahnl2 000). Cri da et al.l d2006t) have generalized the gap 
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opening criterion to simultaneously take into account the mass of 
the planet, the scale height of the disk and the disk viscosity: 
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where H is disk scale height, Rh is the Hill radius, R ey is the ef- 
fective turbulent Reynolds number and q is the mass ratio of planet 
mass (M p ) to stellar mass (M*): 
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Here r p is the planet's semi-major axis, v is the disk viscosity 
and fl p is the angular rotation rate of an object in a circular or- 
bit at r v . This cri terion combines two earlier conditions; the vis- 
cous condition by Brvden e t al. ( 1999), and the tidal condition by 
llin and Papaloizoul i 1993h ! 

Gap opening has been studied exte n sively numerically 
with simulations (e.g. iBrvden etal] 1 19991 ; ICrida et alj 120061 ; 
de Val-Borro et ai] |2006h . However previous work has not pre- 
dicted what mass planets are capable of opening gaps in circumstel- 
lar disks with different accretion rates or structure. The gap opening 
criterion depends on the temperature profiles of the disk through 
the dependence of the viscosity, u, on the sound speed, c s , and ver- 
tical s cale heigh t, H. Here w e adopt an a form for the viscosity 
dShakura and Svu nvaev 1973), 



ac s H. 



(2) 



Since accretion disc viscosity is poorly understood, v and hence 
R ey are not precisely defined in this work; this problem is common 
to all work on accretion discs. 

The disk mid-plane temperature depends on the source 
of heating. We consider viscous heating and heating from ab- 
sorption of radiation from the central star. Passive or non- 
accreting circumstellar disk models that include the effect of 
radiation from th e central star absorbe d onto t he d is k hav e 
been studied by IChiang and Goldreichl dl997l) ; [Belli dl999l) ; 
iDullemond and Dominikl d2004l) . Because of flaring, the irradia- 
tion of the disk star dominates the disk heating at large radii 
when both viscous heating and that from irradiation are considere d 
dCalvet etal.lll99ll : iD'Alessio et ai]|200ll ; iGaraud and Lm]|2007h . 
Passive self-shad owed disks can be considerably colder and thin- 
ner at large radii (Dullemond and Dominik 2004). 

In Section [2] we describe simple disk models and our proce- 
dure for estimating disk temperatures, scale heights. From these, 
we can derive the minimum gap-opening planet mass as a func- 
tion of radius when heating is due to viscous dissipation and when 
heating is due to stellar irradiation. Both cases are considered sep- 
arately so we can be sufficiently flexible to discuss the case of a 
self-shadowed disc. A disc that becomes self-shadowed will have 
significantly reduced heating due to stellar irradiation but heating 
from viscous dissipation would still be present. Section [3] presents 
our results, and a discussion follows in Section|4] 



2 GAP-OPENING IN ACCRETION DISKS 

In th is section we apply t he gap opening criterio n l Crida et all 
\200dt) to the Q-disk model of Shakura and S vunvaevl d 19731) includ- 
ing the effect of stellar irradiation. 

The accretion rate M of a steady thin disk can be calculated 
using mass and angular momentum conservation 



where E is the disc surface mass density. As the disk viscosity de- 
pends on the disk temperature we must consider sources of heat 
to compute it. We use the a prescription to compute the viscosity 
(equation^ but compute the sound speed and disk scale height us- 
ing the disk midplane temperature and the relation for hydrostatic 
equilibrium, 



H = c s /n, 



and the sound speed, 
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where fcs is the Boltzmann constant, /i is the mean molecular 
weight and mn is the mass of the hydrogen atom. In our calcu- 
lation, we take the mean molecular weight [i — 2.4 appropriate for 
interstellar gas. 

We first compute the structure of the disk taking into account 
heat generated from viscous dissipation. Then we consider the case 
of heat generated from the radiation absorbed from starlight. 

2.1 Heat generated from viscous dissipation 



Dissipation due to viscosity gives the energy relation, 
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v = ecri„„, 



(6) 



where e is the emissivity, a is the Stephan-Boltzmann constant, and 
T Vl s is the temperature at the disk surface. A vertical average re- 
lates the surface temperature to the midplane temperature, T VjC by 
Tu,c = f T Tv,s where r is the optical depth from the surface to 
the midplane, r = k S/2, and k is the opacity (e.g., section 2 of 
lArmitage etafll200ll) . We note that the relation between T u s and 
T„, c in Equation[6]is only valid for optically thick disks. In an op- 
tically thin region, we set T„. c to T v , 8 . 

We adopt a convenient analytic form of the emissivity and 
opacity laws that is based on the assumption that dust grains gov- 
ern the opacity and emissivity, e = (T/T©) 6 and k = kv(T /7 1 ©) 6 
where kv=1 cm 2 g -1 , Tq is the solar effective surface temperatur e 
and b = 1 dChiang and Goldreichl [l997l ; IGaraud and Linll2007l) . 
This form is expected as the dust temperature over much of the 
disk corresponds to a peak wavelength (for a black body spectrum) 
that is larger than the diameter of the dust particles and the diam- 
eter is of order the peak wavelength of a Solar temperature black 
body (e.g. lBackman et al.ll 19921 ; IChiang and Goldreichll 19971) . Us- 
ing Equations [2j, [3] — [6] and the opacity law, we solve for the mid- 
plane temperature finding 
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The opacity drops to the low opacity regime at an approximate di- 
viding line of r = 8/3 at a radius 
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For the a disk model, using the above equations, we solve for 
all disk variables (temperature, density, viscosity, scale height) as 
functions of parameters a, M, and M» and as a function of radius, 



M = 37tEia 



(3) 
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2.2 Heating due to irradiation from the star 

The disk can be heated not only by disk viscosity but also by the 
stellar radiation. In this section, we consider the latter case only. If 
we regard the central star as a point source, then the surface tem- 
perature depends on the slope of the disk, 

L» ,„ /dlnff 
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jFranketal.ll2002h . where L, is the stellar luminosity, /3 is the 
albedo and Xi. s is disk surface temperature by irradiation. 

Here, we assume that the disk temperature does not strongly 
depend on height, though a more detail ed model would more 
carefully compute the disk structure (e.g. iD'Alessio et al.l 1 19981 ; 
iGaraud and Linl 12007). Indeed, if the disc is very optically thick, 
then the mid plane temperature will always be dominated by viscous 
dissipation iDu bus et alJll999T) . Ignoring the vertical temperature 
gradient, we estimate the disk midplane temperature Tj, 
using Equations[3]— [6]and|9] 
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where the coefficient 
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We restrict our solution for the disk scale height to a self- 
consistent power law form with din H /dlnr equal to a constant. 
We find that the disk scale height 



H 



u{i-p)T Q y 

127TO- 



(GM*)"9r3, (11) 



and that d In H/d In r = 4/3. If we consider an albedo in the range 
< (3 < 0.5 JWood et al.ll2002h . then Equations [Tol and [TT1 imply 
that the disk structure is not strongly sensitive to the albedo. In our 
subsequent estimates, we have adopted an albedo of (3 — 0. 

Interior to a particular radius, the heat from viscous heating 
dominates that from stellar irradiation. Setting the heating rate due 
to accretion equal to that due to irradiation (setting the left hand 
side of Equation[6]to the left hand side of Equation[l0]and solving 
for radius) we estimate a transition radius 
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For r < rtr we expect that viscous heating would dominate. 
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Figure 1. a) Disk radius vs. midplane temperature for disk mass accretion 
rates M = 10 -9 , 10 -8 and 10 -7 M Q /yr for a = 0.01 and stellar mass 
Mit = IMq. The solid lines are for the case of viscous heating only and 
the dotted lines are for the case of irradiated heating only and L„ = Lq. 
A comparison between the dotted and solid lines shows that heating due to 
stellar radiation dominates at large radius. A change in slope in the solid 
lines occurs at the radius where the disc becomes optically thin (r T ). b) 
Similar to a) except the surface density profile is shown, c) Similar to a) 
except the aspect ratio H/r is shown. 



3 RESULTS 

In Figure [T] we show disk variables as a function of radius for dif- 
ferent accretion rates M = 10" 7 , 10" 8 , and 10" 9 M Q /yr but fixed 
q = 0.01 and stellar mass, M* = IMq. Figure[TJi shows the disk 
mid-plane temperature, Figure[Tp shows the surface density profile 
and FigureQJ the vertical scale height for these disks, and as a func- 
tion of radius in each case. In our Figures, we show with solid lines 
the variables calculated for heating primarily from viscous dissipa- 
tion while variables primarily due to the stellar radiation are shown 



with dotted lines. A comparison between the dotted and solid lines 
shows that heating due to stellar radiation dominates at large radius 
whereas heating due to viscous dissipation dominates at small ra- 
dius. In Figure [TJi-c the change in slope in the solid lines occurs 
where the disc becomes optically thin at a radius approximately 
given by r T (Equation [8}. 

We solve for the minimum gap opening planet mass ratio us- 
ing Equation[T] With a variable substitution of y = (go /<?) 3 equa- 
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tion[T]can be written as a cubic equation 

By + y A = 1 

where we have defined 

50 _ 3 / 3h \3 

qo = Rr y > B = 4-[m^) 

The cubic equation has only one real root0 

B 

with 

3 1 fl ~EP 

u = h \ - H . 

2 V 4 27 

For the entire parameter space covered here we find that B < 1 and 
the minimu m gap openin g planet ITltlSS Qmin ~ qo. Although the 
prediction of lCrida et al.l (Equation[T]in this paper) is more general, 
we find that for most disc models the viscous condition of previous 
work suffices to estimate the minimum gap opening planet mass. 
We wo uld not expect t his to be the case in discs containing 'dead 
zones ( Gammi el 1 1996b . where the viscosity is low, but the scale 
height remains large. 

Figure [2] shows the minimum gap-opening planet mass ratio 
for these disk models as a function of radius. In the inner regions, 
where viscous heating dominates, the minimum q required to open 
a gap is not strongly dependent on radius. However, at larger radii, 
the minimum gap opening q begins to rise more steeply with in- 
creasing radius. This happens for both self-shadowed and irradi- 
ated discs. If the disc is self-shadowed, then we would expect the 
minimum planet mass for gap opening always to lie along the solid 
(viscous heating) line. As noted in Section [2] Dubus et al.l Jl999h 
showed that the midplane temperature of a very optically thick irra- 
diated disc would be dominated by viscous dissipation. In this case, 
the minimum planet mass required to open a gap would lie some- 
where between the solid (viscous) and dotted (irradiated) lines. 
However, unless the disc is self-shadowed, the curve must rejoin 
the dotted line around r T , since then the radiation will be able to 

penetrate to the midplane. 

Figure[2]may be compared to figure 7 of Menou and Goodman 

They were particularly interested in Type I migration rates 
on T Tauri discs, and as part of this, estimated the minimum gap 
opening planet in their models. Our model is showing similar be- 
haviour, with the minimum mass required for gap opening rela- 
tively flat out to ~ 10 au, and rising thereafter. This is despite 
the numerous differen ces in the disc and dust models, a nd the gap 
opening criteria used. lJang-Condell and SasseloV ( 2005) also noted 
that the minimum planet mass required to open a gap would be 
lower nearer to the star. 

Because the minimum gap opening planet mass ratio is less 
sensitive to radius in the inner regions, the transition radius (rt r , 
Equationll2t. between viscous heating and heating by stellar radi- 
ation sets a favourable spot for a planet migrating inward via type 
I migration to slow its migration rate. Type II migration, follow- 
ing gap opening, is expected to be slower than type I m igration for 
Earth mass objects, lacking a gap (Ward and Hahn 2000). A planet 
migrating in the outer disk could become able to open a gap as it 
moved inward. The planets that survive would perhaps be the ones 
that were just massive enough to open gaps in the disc at this transi- 
tion radius; planets unable to open gaps would continue to migrate 

1 See http : / /mathworld .wolfram . com/CubicFormula . html 
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Figure 2. The minimum gap-opening planet mass ratio as a function of 
radius (computed using Equation [T} for the disks shown in Figure [T] with 
accretion rates M = 10~ 7 , 10 -8 , and 10 -9 MQ/yr. The transition radius 
is evident where the dotted line intersects the solid lines. This radius is 
where midplane temperature from viscous dissipation is similar to that from 
stellar radiation. Inside this radius the minimum gap-opening planet mass 
is not strongly sensitive to radius. Outside this radius a larger planet mass is 
required to open a gap, as long as the disk is sufficiently flared to be heated 
by starlight. Lower mass planets can open a gap at the larger transition 
radius (slope change of the solid lines) in self-shadowed disks. 



rapidly, and be lost onto the star. If a flared disk does not contain a 
gap, then the minimum gap opening planet mass ratio computed at 
this transition radius, rt r (Equationll2t provides an estimate for its 
minimum mass. 

As the minimum gap opening planet mass is approximately 
equal to qo, we can compute how it scales with radius, stellar mass, 
luminosity, accretion rate and a parameter. For heating from vis- 
cous dissipation and an optically thick disc, the minimum gap open- 
ing planet mass ratio is 

( M*X~^ ( r \to 
\M Q ) VAU/ 

The minimum gap opening planet mass ratio is only weakly depen- 
dent on radius (with exponent 1/10) and most strongly dependent 
on a and M» . For heating from viscous dissipation and an optically 
thin disc, the minimum gap opening planet mass ratio 




The minimum gap opening planet mass ratio is more strongly de- 
pendent on radius when the disk is optically thin than optically 
thick. 
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When heating is due to irradiation by the star, we find 
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The radial exponent here is 2/3, larger than for either the high and 
low opacity viscous heating cases. For the irradiated disk q m in de- 
pends upon the a parameter because the gap opening criterion de- 
pends on the disk viscosity, however q m m does not depend on the 
accretion rate M because the disk temperature is independent of 
M. 

Were the disk to become self-shadowed then the solid lines 
on Figure[2]and subsequent figures are relevant. The minimum gap 
opening planet mass does not significantly change, but the radius 
at which a low mass object can open a gap is further out and set in- 
stead by the location of r T (Equation[8]l. If the disk is self-shadowed 
then the radius at which the disk becomes optically thin, r T , might 
set a favorable spot for an inward migrating core to open a gap. 

We compute the minimum gap opening planet mass ratio for 
an optically thick disc at the transition radius r± r finding 
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This can be used to place a constraint on possible planets residing 
in a gapless disk, 
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where Mj is a Jupiter mass. The above minimum gap opening 
planet mass ratio and mass are appropriate for both flared and self- 
shadowed disks because the minimum gap opening planet mass is 
insensitive to radius within the transition radius. 

In Figure [3] we show the minimum gap-opening planet 
mass ratio for three disks with the same accretion rate, M = 
10" 8 M o /yr, and stellar mass, M, = 1M Q . but for three differ- 
ent a parameters, a — 0.1,0.01 and 0.001. Figure[3]can be di- 
rectly compared to Figure [2] as the central solid and dotted lines 
are the same. We see that lower mass planets can open a gap in 
lower a disks. We note that the minimum gap opening planet mass 
is more strongly sensitive to this poorly constrained parameter, a 
than the other parameters, such as stellar mass, luminosity and ac- 
cretion rate. 

In Figure [4] we show the minimum gap-opening planet 
mass ratio for three disks with the same accretion rate, M = 
lO -8 M0/yr, and a = 0.01 but with different stellar masses, 
M* = 0.5, 1.0 and 2.0Mq. Figure [4] can be directly compared 
to Figures [2] and [3] The luminosity of a star depends on the stellar 
mass, and here we have assumed that L, = Lq{M*/Mq) to take 
this into account. Figure [4] shows that the lower mass stars (rather 
than higher mass ones) require higher planet mass ratios to open a 

gap- 
Figure |5] shows that the minimum planet mass ratio for three 
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Figure 3. Similar to Figure|2]except the minimum gap-opening planet mass 
is shown as a function of radius for a disk with a = 0.001, 0.01, 0.1, with 
accretion rate M = 10 -8 Mg/yr and stellar mass M, = IMq. The solid 
lines are for the case of viscous heating only and the dotted lines are for the 
case of irradiated heating only. 



0.01 



0.001 



1e-04 




Figure 4. Similar to Figuref2]and[3]except the minimum gap-opening planet 
ratio for different stellar masses Af» = 0.5, 1, 2 Mq for a disk with a = 
0.01 and accretion rate M = 10~ 8 Mq /yr. The solid lines are for the case 
of viscous heating only and the dotted lines are for the case of irradiated 
heating only. We have assumed that L* = Lq (M* /Mq) 3 . 



different mass stars with the additional requirement that the ac- 
cretion rate is proportional to the stellar mass M oc M l, as sug- 
gested by observational surveys dMuzerofle et alj[2005h . With this 
assumption we find that the minimum gap opening planet mass ra- 
tio is approximately independent of stellar mass. Assuming that 
stellar luminosity scales with mass, L* = Lq(M/M,) 3 and 

M = 10~ 8 Af Q /yr (|f^) - Equation [l7]becomes 



M v < 0AM, 



a 
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leaving only a as an undetermined parameter. If larger planets are 
formed earlier during epochs of higher accretion around more mas- 
sive stars then we would predict that they might be more likely to 
be located at larger radii. This follows because the transition radius 
is larger for higher mass stars (Equation [12} as is r T (Equation [8} ■ 
The favorable location for an inward migrating planet to mi- 
grate inward for a flared disk assuming L, = Lq (M/M,) 3 and 

M = lO _8 M0/yr (^jj^^j wou ld be the transition radius com- 
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Figure 5. The minimum planet mass ratio for stellar masses M* = 
0.5, 1, 2 Mq for disks with a = 0.01. We have assumed that M = 
10" 8 M /yr (M*/M ) 2 and L, = L (M*/M ) 3 . The solid lines 
are for the case of viscous heating only and the dotted lines are for the case 
of irradiated heating only. 



puted for these assumptions 
(modifying Equationll2t and that for a self-shadowed disk 
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(modifying Equation [8} Though the minimum gap opening planet 
mass is not significantly different in the self-shadowed disk than 
the flared one, lower mass planets can open gaps at larger radius. 

If larger planets are formed earlier during epochs of higher 
accretion around more massive stars then we would predict that 
they might be more likely to be located at larger radii. This follows 
because the transition radius is larger for higher mass stars. 



4 DISCUSSION AND SUMMARY 

In this paper we have explored the relation between the minimum 
gap-opening planet mass and disk parameters, accretion rate, the 
alpha parameter and stellar mass. We have used simple models for 
q discs, with temperature set by heating due to viscous accretion 
and stellar irradiation. 

We find that the minimum gap opening planet mass is not 
strongly sensitive to radius at radii where the disk is heated pri- 
marily by viscous dissipation and is optically thick However, the 
minimum gap opening planet mass increases with radius where the 
disk is heated primarily by stellar radiation. Due to the weak de- 
pendence with radius in an optically thick disk heated by viscous 
dissipation, the minimum mass gap opening planet in any disk is 
best estimated by the viscous criterion (Equation 1 13t. For a self- 
shadowed disk, the minimum gap opening planet mass is similar to 
that for a flared disk, however the radius at which the mass begins 
to increase more rapidly is at larger radius and is approximately 
where the disc becomes optically thick. 

We estimate that a 0.4 Jupiter mass planet is required to open 
a gap in a flared disk with M = lO _8 M /yr and a = 0.01 around 
a IMq star. Lower mass planets can open gaps in disks with lower 
accretion rates and lower a parameters around a lower mass star. 



We estimate that the minimum gap opening planet mass is propor- 
tional to M OA8 a - s M° A2 LZ om (Equation[P7j. This scaling rela- 
tion can be used to place limits on planets residing in gapless disks 
as a function of accretion rate, a parameter, stellar mass and stellar 
luminosity for both flared and self-shadowed disks. 

It is interesting to speculate on scenarios for planet migration. 
A planet that is not sufficiently massive to open a gap but con- 
tinues to accrete an d migrate could open a gap if it accretes mass 
sufficiently rapidly ( Thommes et al. 2 0071) or because it reaches a 
location in the disk where migration ceases jMassetetai]|2006l : 
ICrida and Morbidellil 120071) . Consider a planet which has formed 
far out in the disk, and is migrating inwards in the linear (Type I) 
regime. If it reaches the transition radius (rv) before achieving the 
mass required for gap opening at that point (Equation 1 16t. then it 
will continue its rapid inward migration, and will probably be lost 
on to the star. The weak dependence of the gap opening mass in the 
inner disk with radius means that it is unlikely that the planet will 
be able to accrete fast enough to open a gap within r tr - The transi- 
tion radius represents a 'last chance' for a planet to open a gap and 
slow its migration. If a disk is observed without a gap, Equationll6l 
provides a bound on the most massive object in the inner portions 
of the disk. 

The different models do predict different locatio ns for planets 
to slow their migration. In the planet trap proposed by Mass et et al.l 
(2006) the truncation of the disk near the star sets the location of 
the planet trap. As accretion would not be detected after a planet 
formed, the scenario suggested here would relate the mass of the 
first planet formed to the minimum detectable disk accretion rate. 
If planet formation cuts off disk accretion and planets stop migra- 
tion at the transition radius then we would predict a relation be- 
tween clearing size and stellar mass. Equation [19] which relates 
the transition radius to the stellar mass assuming L» oc M* and 
M oc M* suggests that the radius of the clearing following the 
first gap opening planet would grow with the stellar mass. We note 
that the radii, r tr and r T , estimated here are of order 1 AU at an ac- 
cretion rate of lO _8 M /yr and so smaller than hole clearing radii 
measured for objects such as CoKuTau/4. The transition radius is 
only as large as CoKuTau/4's (~ 10AU) at an accretion rate of 
order M ~ lO~ 7 Af /yr. 

The disk models considered here are simplistic and do not 
in clude the radiative tran s fer of more sophistica ted models such 
as D'Alessioet ai] J200lh : iGaraud and Linl d2007l) which could be 
used to improve upon the accuracy of our constraints on the mini- 
mum gap-opening planet mass. Better observational constraints on 
the a viscosity parameter would also improve these estimates. We 
have explored disks with opacity with differ ent temperature law s 
for gas disks us ing the opacities described bv lBell and Linl ( 119941) ; 
lBell etal.1 dl997t) The shapes of the curves differ quantitatively 
but not qualitatively from those shown here. In each temperature 
regime the disk temperature has a different slope but the overall 
shape of the temperature, aspect ratio and density curves are simi- 
lar to and within a factor of a few of those shown here. We have not 
shown these curves here as the dust opacity is expected to dominate 
that of the gas, and the simple model explored here can be solved 
more easily to produce our scaling relations. 

Here we have only considered models of steady state disks 
with constant accretion rate M. A starved disk that has not enough 
mass at large radius to maintain its accretion rate would have lower 
disk density at larger radius than predicted here. If the disk were 
flared, the disk temperature would be set by stellar irradiation and 
the minimum gap-opening planet mass estimate would not vary 
from what is predicted here. However were the disk to become self- 
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shadowed the minimum gap opening planet mass would be lower 
at than that predicted with a constant M (this follows because q m in 
decreases with decreasing M; Equations! 13land [74l (. 

Assuming that the accretion processes within the disc can be 
fully parameterised by a is questionable. The source of accretion 
disc viscosity is poorly understood, but the most likely candidate 
is the magneto-rotational instability (MRI). IWinters et alj ( 2003) 
performed magneto-hydrodynamic (MHD) calculations of a planet 
embedded in a disc, and found that the gap structure was markedly 
different to that obtained from a calculation which included a phys- 
ical viscosity. We have also negl ected the possibility of a 'dead 
zone' in the disc dGammiei ri996). where the high midplane den- 
sities shut down the MRI. The surface layers continue to be ionised 
by cosmic rays, and can still accrete. Within the dead zone, the ef- 
fective viscosity is expected to be dramatically lower, and hence 
much lower mass planets might be able to open gaps. In this case, 
the first (tidal) term of Equation [Tj would become more significant 
in determining whether a planet can open a gap. This is in contrast 
to the results we have presented here, which are dominated by the 
second (viscous) term. iMatsumura and Pudritzi 12005. 200(f) con- 
sidered gap opening in a disc with a dead zone. Figure2ofthe l2005l 
paper shows their ex pected minimum mass ga p opening planets, as 
a function of radius. iMatsumura and Pudritz found that the mini- 
mum mass planet required to open a gap was generally increasing 
with radius (unlike the nearly constant value in the inner regions 
we calculated above), with a step at the boundary of the dead zone. 
There are numerous differences between their work and ours; not 
only are the disc models differe nt, but a l so the gap opening cri- 
terion. A criterion calculated by Rafikovl J2002h is used, which is 
based on a calculation of spiral shock dissipation in circumstellar 
discs. The structures formed by a planet which satisfies the vis- 
cous, but not the tidal, condition for gap opening (that is, the first 
term of Equation [T] is greater than 1, but the second term is much 
less than 1) have not been extensively studied numerically, and we 
intend to investigate them further in future work. It should be noted 
that 'dead zones' might not be dead as first thought. Although the 
MRI cannot operate (due to the low ionisation levels), turbulence 
can d i ffuse downwards from the surface layers jFleming and Stonel 
l2003l ; lQishietai]|2007h . stirring material and driving accretion. 

In this paper, we have calculated the minimum mass planet 
capable of opening a gap in a steady state accretion disc. In the 
inner regions, the minimum mass required is very weakly depen- 
dent on radius. Further out, the minimum mass grows with radius. 
The transition between the two regimes (located where heating by 
stellar radiation dominates, or where a self-shadowed disc becomes 
optically thin) represents a 'last chance' for a rapidly migrating em- 
bedded planet to open a gap. If a planet is unable to do this, it is 
likely to be lost onto the star. 
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